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ABSTRACT 

We present multi-molecular line maps obtained with the Mopra Telescope towards the 
southern giant molecular cloud (GMC) complex G333, associated with the Hll region 
RCW 106. We have characterised the GMC by decomposing the 3D data cubes with 
GAUSSCLUMPS, and investigated spatial correlations among different molecules with 
principal component analysis (PGA). We find no correlation between clump size and 
line width, but a strong correlation between emission luminosity and line width. PGA 
classifies molecules into high and low density tracers, and reveals that HCO+ and 
N2H+ are anti-correlated. 

Key words: stars: formation - ISM: clouds - ISM: molecules - ISM: structure - radio 
lines: ISM 



1 INTRODUCTION 



This is the third paper (paper I by Bains et al. 



2006 



and 



paper II by I Wong et al.|2008t in a series of multi-molecular 



lines observations of the giant molecular cloud (GMC) 
associated with RCW 106, G333. The G333 giant molecular 
cloud, a southern massive star forming region, spans 

1?2 X 0?6 on the sky roughly centred on / - 333°, b 0?5 

(aj2ooo = 16^21^, 5j2ooo = -50^30"^), at a distance of 



3.6 kpc (Lockman 1979). The aim of the multi-molecular 



line study was to investigate the relationship between the 
dynamics of the interstellar medium and star formation. 
Further investigation and analysis planned include power 
spectra to study the role of turbulence in the GMC, and 
searching for evidence of triggered star formation. 



The G333 GMC has previously been studied in other 
wavelengths, for example in the far-infrared by |Karnik 



et al. (2001) and the 1.2-mm dust continuum by Mookerjea 



et al. (2004). A search for water masers was carried out 



by |Breen et al.| (|2007|. There are also numerous works 

(e.g. 



on specific regions in this GMC, such as RCW 106 
Rodgers, Campbell Whiteoak||T960l [Russeil et a"L||2005 ) 



and G333.6— 0.2, one of the brightest southern compact 
Hii regions (e.g. Goss k Shaver 1970 , Becklin et al.||1973[ 



Store y~eFlal|1989| |Colgan et al.||1993| [Fujiyoshi et al.| 



,1998| |2001| |2005| and |2006|). However, to date, there has 



been no systematic, 3-mm multi-molecular line mapping 
of this GMC. This study demonstrates the full capabil- 
ity of the Mopra Telescope'^ new digital filter bank. 



E-mail: nlo@phys.unsw.edu.au 
t Bolton Fellow, ATNF 



^ The Australia Telescope Mopra telescope is part of the Aus- 
tralia Telescope, which is funded by the Commonwealth of Aus- 
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Table 1. List of observed molecular transitions mapped in G333. The columns are: (1) molecule; (2) transition; (3) rest freq uency (|Lovas,| 
Johnson Snyder|^19 79 ); (4) upper energy levels from the Cologne Database for Molecular Spectroscopy (CDMS, M tiller et al.||200Tf 
2QQ5| ); (5) calculated critical density: unless specified, ricrit = ^nz/(<^('^)'^} where A^i is the Einstein A coefficient obtained from CDMS, 
{(j{v)) = lO"-*^^ cm-^ is the collision cross section and v is the velocity of collision particles assumed to be 1 km s"-*^; (6) observation 
season; (7) backend used, where MPCOR stands for Mopra Correlator, which was decommissioned in 2005 and replaced with a new 
digital filter bank, the UNSW-Mopra Spectrometer (UNSW-MOPS); (8) notes, where [hf] denotes hyperfine structures, with the rest 
frequency quoted being that of the strongest hyperfine component, and [f] denotes molecules that have detectable emission at a few 
places only. Their emission maps are not presented in this work but will appear in subsequent papers on specific regions in the G333 
cloud. 



Molecule 


Transition 


Rest frequency 


Eu/k 


T^crit 


Observation season 


Backend 


Notes 






(GHz) 


(K) 


(cm-3) 








13CO 


1-0 


110201.353 


5.29 


3 X 103 


2004 


MPCOR 




ci^o 


1-0 


109782.173 


5.27 


3 X 10^ 


Jul, 2005 


MPCOR 




CsH 


1-0 


87316.925 


4.19 


2 X 10^ 


Jul, 2006 


MOPS 


[hf] 


CH3OH 


2(0,2)-l(0,l)A+ 


96741.377 


6.97 


7 X 10^ 


Sep, 2006 


MOPS 


[t] 


cs 


2-1 


97980.953 


7.05 


2 X 10^ 


Sep, 2006 


MOPS 




C34S 


2-1 


96412.961 


6.94 


2 X 10^ 


Sep, 2006 


MOPS 


[t] 


HCN 


1-0 


88631.847 


4.25 


2 X 10^ 


Jul, 2006 


MOPS 


[hf] 




1-0 


86340.167 


4.14 


2 X 10^ 


Jul, 2006 


MOPS 


[hf,t] 


HCCCN 


10-9 


90978.989 


24.01 


6 X 10^ 


Sep, 2006 


MOPS 


[t] 


HCCCN 


11 - 10 


100076.385 


28.82 


8 X 10^ 


Sep, 2006 


MOPS 


[t] 


HCO+ 


1-0 


89188.526 


4.28 


4 X 10^ 


Jul, 2006 


MOPS 




H13CO+ 


1-0 


86754.330 


4.16 


4 X 10^ 


Jul, 2006 


MOPS 


[t] 


HNC 


1-0 


90663.572 


4.35 


3 X 10^ 


Jul, 2006 


MOPS 


[hf] 


N2H+ 


1-0 


93173.480 


4.47 


4 X 10^ 


Sep, 2006 


MOPS 


[hf] 


SiO 


2-liy = 


86847.010 


6.25 


3 X 10^ 


Jul, 2006 


MOPS 


[t] 


SO 


3(2) -2(1) 


99299.905 


9.23 


1 X 10^ 


Sep, 2006 


MOPS 


[t] 




4(5) -4(4) 


100029.565 


38.58 


3 X lO'Q 


Sep, 2006 


MOPS 


[t] 



Ungerechts et al. (19971 



the UNSW-Mopra Spectrometer (UNSW-MOPS). Other 
similar surveys utilising the UNSW-MOPS are the H2O 
Southern Galactic Plane Survey (HOPS) by Walsh et al. 
(2008) and the Central Molecular Zone of the Galaxy by 



Jones et al. (2008) 



We present in this paper results of molecular line 
mapping carried out during July to November in 2006, con- 
sisting of molecular lines from 83 to 101 GHz. For a detailed 
analysis of ^^CO and C^^O data please refer to previous 
work by 'Bains e t al.| ( [2006 , hereafter BWC2006) and | Wong 
et al.| f2QQ8, hereafter WLB2008). In Section |2] we describe 



the observing technique, and the new UNSW-MOPS digital 
filter bank. In Section [3] we present the results, and examine 
the velocity and spatial distribution of the molecular 
emission from the different molecules. In Section [4] and 
[5] we present two different approaches to characterising 
emission distribution in the GMC: Clump finding with 
GAUSSCLUMPS and principal component analysis (PGA). In 
Section [6] we discuss results, including specific regions of 
the GMG. Finally, we summarise our findings in Section jT] 
Data presented here will be made available in the near 
future, please contact the authors for details. 



tralia for operation as a National Facility managed by CSIRO 
(ATNF, http: /7www.atnf.cs iro.au| ). 



2 OBSERVATIONS 

The data were collected with the 22 metre Mopra Telescope, 
which is a centimetre- and millimetre-wavelength antenna 
having a full width to half-maximum (FWHM) beam size 
of -36''at 100-GHz (iLadd et al.||2005|). The observations 



were carried out with the narrow band mode of the new 
UNSW-Mopra Spectrometer (UNSW-MOPS) digital fil- 
terbank back-end, and a Monolithic Microwave Integrated 
Circuit (MMIG) 77 to 116 GHz receiver. The observing 
parameters, including rest frequencies and dates are listed 
in Table [l] UNSW-MOPS has a 8-GHz bandwidth with 
four overlapping 2.2-GHz subbands, each subband having 
four dual-polarisation 137.5-MHz-wide windows giving a 
total of sixteen dual-polarisation windows. Each window 
has 4096 channels providing a velocity resolution of ~0.1 
km s~^ per channel at 100 GHz. The velocities presented in 
this work are with respect to the kinematic local standard 
of rest (LSR), with t'LSR = — 50 km s""*^ being the systemic 
velocity of the G333 complex. 

The brightness temperature Tb is related to the 
antenna temperature by Tb = TX/?7i/, where rfu is 
the frequency dependent beam efficiency. According to 



Ladd et al. (2005) the main beam efficiency at 86 GHz is 
7786 GHz = 0.49, and at 110 GHz is ?7iioghz = 0.44. The 
beam efficiencies are ?786GHz = 0.65 and ?786GHz = 0.56 
at 86 and 110 GHz respectively. The results presented in 
this paper are in terms of antenna temperature TJ unless 
otherwise specified. 

We followed a similar mapping procedure to our 
previous work (^^CO; BWC2006). The -1 square degree 
region was divided into 300 x 300 square arcsec fields, 
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with pointing centres separated by 285 arcsec in right 
ascension and dechnation, ahowing a 15-arcsec overlap 
between adjacent fields. Each field was mapped with two 
passes, the first pass scanning in right ascension and the 
second pass scanning in declination. The observing mode 
was 'on-the-fly'(OTF) raster scanning, at a scan rate of 
3.5 arcsec s~^ and averaging data over a 2-s cycle time. 
We used the ^^CO map as a guide, mapping the brightest 
^^CO fields first, then extending these to cover the majority 
of the emission in most detected molecules, resulting in 
maps covering ^0.7 square degrees. We used two frequency 
settings centred at 87 and 97 GHz, with each setting 
covering an 8-GHz bandwidth, 4-GHz on either side of the 
central frequency. 

We have used the SiO masers IRS VI 540 and AH 
Scorpii for pointing calibration. Pointing calibration was 
performed at the beginning of each map, therefore approx- 
imately every 50 minutes. Estimated pointing errors are 
within 10 arcseconds. 

At the end of each day's observation, a 5 minute posi- 
tion switched observation of M17SW was carried out with 
the same frequency settings as for the maps, to monitor 
the performance of the system. Based on these observations 
the deviation of during the whole observing season is 
within 10 per cent. 

The data were reduced with the livedata and 
GRIDZILLA packages available from the ATNF, written by 
Mark Calabrettc[j livedata performs a bandpass calibra- 
tion using the preceding reference (off-source) scan, then 
fits the spectral baseline with a first degree polynomial. Due 
to low sensitivity near the start and end of the window, the 
first and last 750 channels were discarded before baseline 
fitting. GRIDZILLA grids the data according to user specified 
weighting and beam parameter inputs. In this work, our 
data was weighted by the relevant system temperature 
(Tsj^s) and Gaussian smoothed during data reduction stage. 
The data cubes were gridded with three pixels across each 
beam width, therefore each pixel is 15 x 15 arcsec^ for 
data with 87-GHz frequency setting, and ~ 12 x 12 arcsec^ 
for the 97-GHz setting. Note that all data cubes presented 
in this work are regridded to 15 x 15 arcsec^ pixel size for 
ease of comparison, and the velocity resolution is left as 
- 0.1 km s"\ 



3 RESULTS 

3.1 Spatial distribution of the molecules 

Different molecules may trace different physical conditions 
within the region, such as density, chemistry, temperature, 
evolutionary stage, and so on. Our selection of molecules 
covers those which may trace differences in density, temper- 
ature and dynamic features such as outflow, allowing us to 
study environmental varieties within this giant molecular 
cloud. Since some molecules (e.g. CH3OH, HCCCN) have 
detectable emission at a few places only, their emission 
maps and velocity profiles are not presented here, but in 
subsequent works focusing on specific regions. 

^ http:/ /www. atnf.csiro.au/computing/software/ 



Table 2. A list of the emission range, peak brightness and Icr 
level of the integrated emission maps shown in Figures [l]and[2] 



Molecule 


Emission 


range 


1 a 


Peak brightness 




(km s~ 




(K km s~^) 


(K km s~^) 


13CO 


-70 to 


-40 


0.3 


95 




-70 to 


-40 


0.2 


19 


CS 


-70 to 


-40 


0.2 


30 


HCO+ 


-70 to 


-40 


0.1 


21 


HCN 


-80 to 


-30 


0.1 


21 


HNC 


-70 to 


-40 


0.1 


18 


N2H+ 


-80 to 


-30 


0.2 


11 


C2H 


-60 to 


-40 


0.2 


5 



Shown in Figure [T] and [2] are the integrated intensity 
(zeroth moment) maps of the observed molecular lines 
(contours) overlaid on the Spitzer Galactic Legacy Infrared 
Midplane Survey Extraordinaire (GLIMPSE; |Benjamiri| 



et al. 2003) 8.0-/jm image (grey scale). The integrated 
velocity range is -70 to -40 km s"^ for ^^CO (BWC2006), 
C^^O (WLB2008), CS, HCO+ and HNC, and -60 to -40 
km s~^ for C2H. A wider velocity range was chosen for 
molecules with hyperfine components: —80 to —30 km s~^ 
for HCN and N2H^, so as to include emission from all 
components. The choice of velocity range reflects emission 
from the GMC only. Emission outside this range comes 
from spiral arms and intermittent clouds. The maps were 
Gaussian smoothed with an FWHM of 45 arcsec. The labels 
in the ^^CO total intensity map (Figure [TJi) indicate some 
of the molecular regions in proximity to bright Hll regions 
and IRAS sources, as marked by rectangles. The starting 
contour levels are based on the lowest detectable emission 
above noise, hence outlining the emission structure of the 
cloud. Note that contours near the edges of the maps are 
noisy; this is because the edges were scanned with one 
pass only and thus have lower sensitivity. Details of the 
integrated emission maps are summarised in Table [2] 

The CS total intensity map (Figure [l]:) shows a more 
confined distribution than the ^^CO emission. This is 
as expected from a dense gas tracer such as CS, with a 
critical density of ~ 10^ cm~^. However CS has a similar 
distribution to another CO isotope presented here — C^^O. 
HCO+ (Figure [ijl), a common ionic species found in 
molecular clouds also presents a similar distribution to CS, 
as do HCN (Figure [2ji) and HNC (Figure [2})). 

A comparison of HCN and HNC total intensity maps 
(Figures [2^ and b), shows that they both have similar 
distributions in general. This is surprising because of 
the strong chemical differences between HCN and HNC, 
meaning that its abundance ratio can depend heavily on, 
for example, temperature (e.g. Schilke et al. 1992 Hirota 



et al. 1998). We do note there are local variations where 
temperature changes, such as the massive dense cold core 
where HCN and HNC have different spatial distribution 
dLo et al.|2007 i). 

The total emission map of N2H+ (Figure ^) shows 
that this molecule has a more compact and clumpier dis- 
tribution than other molecules. It also shows regions with 
intense N2H^ emission which show diffuse emission in other 
molecules; here we highlight two such regions on the map 
(indicated by arrows). The source G333. 125—0.56 (bottom 
arrow) is a massive dense cold core, with temperatures of 
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Figure 1. Integrated emission maps (contours) for ^^CO {J = 1- 0), C^^O ( J = 1 - 0), CS ( J = 2 - 1) and HCO+ ( J = 1 - 0) overlaid 
on a Spitzer IRAC 8.0-/im image (grey scale). The maps are integrated over a velocity range of -70 to -40 km s'^ and were clipped at 
a 3 a level. 1 cr is 0.3 K km s'^ for I^CO, 0.2 K km s'^ for C^^O and CS, and 0.1 K km s'^ for HCO+. The maps were then Gaussian 
smoothed with an FWHM of 45 arcsec. The contour levels start at 1.5 K km s-^ for ^^CO, 2.0 K km s-^ for C^^O, 1.4 K km s^^ for CS 
and 1.0 K km s~-^ for HCO+. Then each successive contour level is double the value of the previous one. The choice of lowest contour 
level is based on the lowest detectable emission above noise, except for -"^^CO where there is emission detected throughout the map. 
The temperatures are in terms of the antenna temperature, TJ. Previously designated Hll regions and the associated IRAS sources are 
marked with white rectangles and names are shown in the -"^^CO map. The scale bar shown in the ^^CO map indicates 5-pc at a distance 
of 3.6-kpc and the tiny circle in the lower left corner of each panel show the FWHM beam size after smoothing. The velocity range used 
here reflects emission from the GMC only; emission outside this range comes from more distant spiral arms or intervening clouds. 
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Figure 2. Integrated emission maps (contours) for HCN {J =1- 0), HNC {J=l- 0), N2H+ (J =1-0) and C2H {J = 1-0) overlaid 
on Spitzer IRAC 8-/im images (grey scale). The integrated velocity ranges are of -80 to -30 km s-^ for HCN and N2H+, -70 to -40 
km s-i for HNC and -60 to -40 km s-^ for C2H. The maps were clipped at a 3 cr level. 1 cr is 0.1 K km s-^ for HCN and HNC, and 
0.2 K km s~-^ for N2H"^ and C2H. The maps were then Gaussian smoothed with an FWHM of 45 arcsec. The contour levels start at 2.0 
K km s-i for HCN, 1.3 K km s'^ for HNC, 1.6 K km s'^ for N2H+ and 1.2 K km s'^ for C2H, then each successive level is double 
the previous one. The choice of lowest contour levels are based on the lowest detectable emission above noise. The temperatures are in 
terms of and the tiny circle in the lower left corner of each panel show the FWHM beam size after smoothing. Previously designated 
Hll regions and the associated IRAS sources are marked with white rectangles (see Figure [T] for their naming). The arrows indicate two 
N2H"^ emission peaks, which are not prominent in other emission maps. Note that the differences in velocity range between molecules 
are due to hyperfine splitting such that some molecules have a wider velocity distribution. The velocity range used here reflects emission 
from the GMC only; emission outside this range is from more distant spiral arms or intervening clouds. 



6 A^. Lo et al 



13 to 19 K derived from the NH3 (J,K) = (1,1) and (2,2) 
inversion line and from the spectral energy distribution 
(SED). It has broad thermal SiO emission, and is believed 
to be a deeply embedded young stellar object in an early 

In the northern 



2007) 



stage of star formation (Lo et al, 
part of the map, the N2H^ emission indicated by the top 
arrow is similarly interesting. Here, bright N2H^ emission 
is detected where other molecules show diffuse weak 
emission. A comparison with the GLIMPSE 8.0-/im image 
clearly shows that N2H^ aligns well wit h the infrared dark 

While N2H+ 



6.3 



filaments, as discussed later in Section 
seems to be correlated with infrared dark filaments, C2H 
does the opposite; it seems to correlate with GLIMPSE 
8.0-/xm emission. 

Among the eight molecules we present, C2H emission 
is the weakest. Its spatial distribution is more extended 
than the quiescent gas tracer N2H^, but not as extended 
as other molecules. 



3.2 Velocity distribution of the molecules 

Spatially averaged spectral emission profiles of the molecules 
presented in Section |3.1| are shown in Figure [3| Note that 
the apparent lower signal-to-noise of the N2ll^ and C2H 
spectra is due to these two molecules being less extended 
spatially. To illustrate this, spectra of specific regions 
(IRAS16172-5028, IRAS16164-5046 and G333.125-0.56) 
are shown in Figure [4] showing that the signal-to-noise 
level is comparable among the molecular lines. From the 
spectral profiles in Figure [3] we can see that the molecular 
emission peaks at ^ —50 km s~^, the systemic velocity of 
the G333 cloud, and spans a velocity range of '-^ —70 to 

40 km s-^ ^^CO, C^^O and CS have similar velocity 

structure, such as the emission 'shoulders' at ~ —40 and 
- -55 km s"\ The - -55 km s"^ emission feature also 
appears in the HCO+ and HNC profiles. HCN and N2H+ 
have the widest velocity range, but this is due to hyperfine 
splitting rather than actual differences in distribution. The 
^^CO velocity feature centred at —70 km s~^ comes from a 
different cloud along line of sight, as discussed in BWC2006. 

From the ^^CO data, BWC2006 noted a linear velocity 
gradient of ~ 0.2 km s~^ pc~^ across the GMC (1 arcmin 
is 1 pc for this GMC at distance of 3.6 kpc), which is five 
times larger than the velocity gradient due to Galactic 
rotation (~ 0.04 km s~^ pc~^ for an angular separation 
of ^65 arcmin at 6 = 333°, see BWC2006). To investigate 
whether molecules other than ^^CO are involved in this 
bulk gas motion, we plotted the position-velocity (pv) map 
of CS, HCO+, HNC and C2H across the GMC as shown 
in Figure [5] The pv slice is positioned such that it cuts 
through the major emission ridge of the GMC, as shown 
in the top panel inset. The centre position (zero angular 
offset) is roughly at aj2ooo = 16^21^, ^j2ooo = -50^30"^ 
near IRAS16172— 5028, with positive offset corresponding 
to the southern part of the GMC. 

Evident in the CS pv map are a number of emission 
clumps; those near the GMC centre (centred on 0' offset) 
are surrounded by diffuse gas. The line widths of these 
emission clumps are 4 to 14 km s~^. Similar to ^^CO, 
the velocity gradient is also noticeable in CS, HCO^, HNC 
and C2H, ^ 10 km s~^ over roughly 50 arcmin {r^ 0.2 km 



s"^ pc"^), comparable to that of the ^^CO. The HCO+ 
pv map shows similar overall structure to CS in general, 
however the centre velocity of the strong emission clump at 
—30 arcmin (G333.6-0.2, as labelled in Figure [sf is -50 
km s~^ for HCO^, compared to —48 km s~^ for CS, with 
HCO^ emission being more extended towards the negative 
velocities (- -60 km s"^) than CS. Where the CS emission 
clump peaks at —48 km s~^, the HCO^ emission clump is 
less intense and appears to host two components. In the 
HNC emission, the structure of this feature is consistent 
with CS. In the 3-D data cube, the HCO^ spectrum 
shows a deep self- absorption feature at this position, when 
compared with the optically thin H^^CO^, explaining the 
absence of HCO^ emission. In general, the HNC and C2H 
pv maps show a similar velocity structure to CS, except 
there is no detectable C2H at ~ 28 arcmin at —58 km s""*^. 

In the next two sections, we will present two different 
approaches in characterising the distribution of the emis- 
sion in the GMC, clump finding with GAUSSCLUMPS and 
principal component analysis. 



4 CLUMP ANALYSIS 

One of the aims of this multi-molecular line mapping is 
to examine how different molecules correlate with each 
other in the GMC. To characterise the distribution of 
molecules, we have used the Gaussian clump decomposition 
algorithm G AUSSCLUMPS ( Stutzki Guesten]p^990 : Kramer] 
et al. 1998) to decompose the observed three-dimensional 



data cubes into individual 'clumps' of emission. There has 
been some discussion in recent years about the best type 
of automatic clump decomposition algorithm to use in 
giant molecular clouds (see Sheth et al. 2008 for a good 
discussion on this subject). GAUSSC LUMPS and clumpfind 
(Williams, de Geus & Blitz 1994) are algorithms which 



work by identifying peak pixels within a data cube or 
image, and then identifying other nearby pixels which 
may be part of the same clump, clumpfind stops adding 
pixels to a clump when it reaches pixels below a user 
specified contour level (see BWC2006 for a more lengthy 
discussion of this algorithm). GAUSSCLUMPS decomposes 
the data cube by iteratively fitting the peak position as a 
Gaussian distribution, then subtracting the fitted clumps 
and fitting the residual map. Both these methods have 
the disadvantage that they tend to break the emission 
into many small clumps close to the resolution element of 
the data set, particularly with a low signal-to- noise data 

different 
which 



set 



Rosolowsky & Leroy (2006) have described 



method, implemented by the algorithm cprops, 
uses a combination of moments and principal component 
analysis to identify clumps in data cubes, which they find 
is more robust against the effects of resolution and noise. 
However, of the three methods, we find GAUSSCLUMPS 
provided the best set of physically plausible clumps for the 
molecular transitions presented in this paper, which trace 
dense confined gas, with CPROPS finding very few clumps 
in each data cube. This is consistent with the finding of 



Rosolowsky & Leroy ( 2006 ) that GAUSSCLUMPS is the most 



effective algorithm for separating tight blends of clouds. 
Consequently, we have chosen to use the GAUSSCLUMPS 
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Velocity (km s""!) 



Figure 3. Spatially averaged emission profiles of the whole GMC. The temperatures are in terms of T^. Note that the lower signal-to- 
noise ratio of the N2H+ and C2H spectra is due to these two molecules being less extended spatially; the 1 cr level is comparable to other 
more extended molecules. 



algorithm, but have also undertaken a careful check of 
the clumps produced to confirm that they are physically 
plausible. 

The molecules selected for such clump finding are 
CS, HCO+, HNC and C2H. The selection is based on 
CS and HCO^ emission not having hyperfine splitting, 
HNC hyperfine splitting not being resolvable, and the C2H 
hyperfine component separation being so wide that only 
one component falls into the velocity range. The data 
cubes were smoothed with a hanning window width of 
five channels and binning window of two channels, giving 
a velocity resolution of ~ 0.3 km s~^, to improve the 



signal-to-noise ratio. The rms levels after binning are 0.070, 
0.065, 0.06 and 0.06 K (in antenna temperature, T^) per 
0.3 km s-^ for CS, HCO+, HNC and C2H respectively. We 
have set the intensity threshold at the 5cr level, therefore 
clumps with peak temperature below 5cr are discarded 
for CS, HCO+ and HNC. For C2H, due to its compara- 
tively weak emission, we only discarded clumps with peak 
temperature below the 3cr level. For all molecules the 
spectral line emission towards the clumps was examined to 
confirm that the emission is real. We also excluded clumps 
that have angular sizes smaller than 1.5 beam- widths to 
reduce the number of false detections. Clumps that are 
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IRAS1 61 64-5046 



IRAS1 61 72-5028 



G333. 125-0.56 




-50 

Velocity (km/s) 



-60 -50 

Velocity (km/s) 



70 -60 -50 -40 
Velocity (km/s) 



Figure 4. Sample spectra towards two of the IRAS sources (IRAS16172-5028 and IRAS16164-5046) and the cold core G333.125-0.56 
in G333, showing the general quality of the data, and the consistent noise level within the data set. Line temperatures are given on a 



found on the edges of the maps are also discarded, due to 
the lower sensitivity as the map edges contain one pass only. 



4.1 Clump properties 

GAUSSCLUMPS finds 129, 186, 128 and 78 clumps in the 
CS, HCO^, HNC and C2H data cubes, respectively. 
GAUSSCLUMPS reports centre positions and velocities, line 
widths, angular sizes, peak temperatures and internal 
velocity gradients of the identified clumps. We have listed 
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Figure 5. Position- velocity diagrams of CS, HCO+, HNC and C2H centred at approximately q;j2ooo = 16^21"^, Sj2qqq — —50^30^. 
The inset in the top panel shows the position-velocity cut presented in the diagrams. The contour levels all start at the 3cr level, 0.3 K, 
increasing to 3.0 K with increments of 0.5 K. The temperatures are in terms of T^. Positive angular offset corresponds to the southern 
region of the GMC. The arrow points toward the Hll region G333.6— 0.2. 
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Table 3. The 20 brightest CS clumps selected according to decreasing peak temperature with line intensities above the Scr level of the 
Hanning smoothed data cube. Also in Tables ^ [s] and [S] The listed properties are centre position (RA and Dec), centre velocity {v in 
km s"-*^), full width to half maximum (FWHM) line width (AV in km s-^), FWHM of the two principal axes (Dx and Dy in arcmin), 
peak temperature (T^ = T^/'q^ in K), internal velocity gradient across the clump {dv/dr in km s"-*^ arcmin" -'^), calculated luminosity (L 
in K km s"-*^ pc^), L = (d [pc])^( i3Q^3gQQ )'^(rxry) f Tij6v (see text for details) and associated clumps in other molecules. See the online 
version for a complete list of 129 clumps. 

# RA Dec V AV Dy Peak Tb dv/dr L Associated clumps in other molecules 



CSOl 


16: 


:21:02.0 


-50:35:04.9 


-52.2 


2.8 


1.7 


1.4 




1 fi' 
±o, 


■91 -01 


-ou.oo.zz.y 


AQ ^ 
-4ty .0 


9 ^ 


±.0 


±.0 




iO, 


,99.07 8 


-ou.uD.zz.y 


/I 7 8 


A 1 


9 n 
z.u 


i.o 




16: 


:21:02.0 


-5U:o4:5o.y 


-55.3 


3.1 


1.0 


1.2 


CS05 


16: 


:20:10.2 


-50:53:16.9 


-57.4 


3.4 


1.5 


1.6 


CS06 


16: 


:21:27.1 


-50:24:52.9 


-50.4 


2.6 


1.9 


2.4 


CS07 


16: 


:21:30.3 


-50:26:46.9 


-52.2 


3.2 


1.4 


1.7 


CS08 


16: 


:21:13.3 


-50:39:46.9 


-55.6 


2.0 


1.3 


1.7 


CS09 


16: 


:20:11.5 


-50:53:22.9 


-54.2 


3.2 


1.1 


0.9 


CSIO 


16: 


:20:26.0 


-50:41:22.9 


-56.0 


2.4 


1.2 


1.4 


CSll 


16: 


:20:46.9 


-50:38:46.9 


-53.3 


2.6 


1.4 


1.2 


CS12 


16: 


:21:37.8 


-50:24:52.9 


-50.8 


3.0 


1.1 


2.4 


CS13 


16: 


:19:37.1 


-51:03:28.9 


-50.6 


3.7 


1.7 


1.2 


CS14 


16: 


:22:07.8 


-50:06:04.9 


-45.3 


2.0 


1.7 


1.4 


CS15 


16: 


:20:09.6 


-50:53:16.9 


-60.1 


2.4 


1.4 


1.4 


CS16 


16: 


:21:30.3 


-50:26:52.9 


-54.8 


2.5 


1.2 


1.6 


CS17 


16: 


:21:16.5 


-50:39:40.9 


-57.6 


2.0 


1.2 


1.6 


CS18 


16: 


:20:38.0 


-50:41:34.9 


-54.8 


2.1 


0.9 


2.0 


CS19 


16: 


:21:18.3 


-50:30:34.9 


-52.2 


2.0 


0.9 


2.2 


CS20 


16: 


:21:13.3 


-50:33:40.9 


-51.5 


2.1 


1.1 


2.0 



the parameters of the brightest 20 clumps of CS, HCO^, 
HNC and C2H emission according to their peak intensities, 
along with the derived luminosity in Table [Sj [4] |5] and |6] 
respectively. A complete list of the clumps is available in 
the online version. The temperatures listed are corrected 
for the beam efficiency, T^ = T^/?7zy, where 7]^ is the 
extended beam efficiency. For the twenty brightest clumps 
of each molecular species, we also listed their associated 
clumps in other molecules in the last column of the 
tables. The association criteria are: the clumps coincide 
within one beam and the line widths overlap. We have 
also plotted the clump positions on the corresponding 
integrated intensity map in Figures [6] and |7| with ellipses 
of size proportional to the angular sizes of the clumps. A 
summary of the clump properties is listed in Table |7| for 
comparison. The centroid velocities and line widths of the 
clumps among the four molecules are similar. The angular 
sizes of the clumps are also comparable among the four 
molecules, but note that the lower limit of angular size is 
just above the threshold value (1.5 beam- width) we have set. 



4.2 Optical depth and column density 

Among the four molecules decomposed by GAUSSCLUMPS, 
CS and HCO^ have lines from rare isotopes, so-called 'iso- 
topologues' (C^^S and H^^CO+) in our data set. Therefore it 
is possible to estimate optical depths for these molecules. As- 
suming identical excitation temperatures and optically thin 
emission from the rare isotopologues, the opacity r relates to 



7 


n A 
U.4 


/I Q 7 


ilLyWpUz, rllNLyUz, LyZrlUi 


^ 9 
<J. z 


n 1 
u. 1 




y\nr)r^r\9 m\rpn9 mvpnf^ 

iiv^wpuz, niNouz, niNvyUtj, vyznui 


A Q 


n A 


/IQ 1 


v^oi^, iiv^wpuo, nv^wpuy, riiNv^ui. 








illN V^iU 


4.4 


0.1 


14.8 


HCOp05, C2H05 


3.8 


0.8 


29.8 


CS09, CS15, HCOpl8, HNCll, 








HNC17, C2H06 


3.6 


0.2 


42.2 


HCOp03, HNC03, C2H08 


3.6 


0.2 


26.7 


CS16, HCOpl6, HCOp20, HNC12, 








C2H02 


2.9 


0.5 


12.4 


CS17, HCOp07, HNC07, C2H07 


2.6 


0.5 


8.1 


CS05, HCOpl8, HNC17 


2.6 


0.3 


9.9 


HCOpM, HNC16, C2H11 


2.5 


0.3 


10.4 


HNC18, C2H04 


2.0 


0.5 


16.1 


C2H15 


2.0 


0.7 


14.6 


C2H12 


2.0 


0.3 


8.5 


CS03, HCOp04, HNCOl, C2H03 


2.0 


0.6 


8.6 


CS05, HNCll 


1.9 


0.6 


9.6 


CS07, HCOpl6, HCOp20, HNC12, 








C2H02 


1.9 


0.1 


7.3 


CS08, HCOp07, HNC07, C2H07 


1.9 


0.3 


6.9 




1.9 


0.4 


7.1 


C2H13 


1.8 


0.4 


8.3 


HNC15 



Table 7. A summary of the clump properties for the CS, HCO+, 
HNC and C2H data cubes, derived from GAUSSCLUMPS (see Sec- 
tion [IJ. 





CS 


HCO+ 


HNC 


C2H 


Number of clumps 129 


186 


128 


78 


Centre velocity (km s~^) 


lowest 


-61.9 


-64.8 


-65.9 


-61.4 


highest 


-44.0 


-44.4 


-44.4 


-43.5 


mean 


-52.0 


-52.6 


-52.6 


-52.4 


cr 


4.1 


3.8 


3.7 


3.5 


(arcmin) 


minimum 


0.9 


0.8 


0.8 


0.9 


maximum 


2.3 


2.0 


2.3 


1.8 


mean 


1.4 


1.4 


1.3 


1.3 


a 


0.2 


0.2 


0.3 


0.2 




Peak Tb 


(K) 






minimum 


0.7 


0.7 


0.7 


0.5 


maximum 


7.6 


5.4 


4.4 


1.8 


mean 


1.3 


1.6 


1.8 


0.7 


cr 


1.0 


0.7 


0.7 


0.3 




dV (km s 








minimum 


0.4 


0.5 


0.4 


0.3 


maximum 


4.1 


3.8 


3.7 


3.2 


mean 


1.7 


1.1 


1.1 


1.5 


cr 


0.7 


0.6 


0.6 


0.6 




Luminosity (K k 


;m s""*^ pc^] 


) 




minimum 


0.05 


0.04 


0.02 


0.03 


maximum 


5.2 


5.0 


4.0 


1.8 


mean 


0.6 


0.4 


0.5 


0.2 


a 


0.9 


0.5 


0.6 


0.3 
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Figure 6. The integrated emission maps (contours and grey scale) of CS and HCO+ (for the contour levels, see Figure ^ overlaid with 
the corresponding GAUSSCLUMPS clump fits of the 3-D data cube. The ellipses show the orientation and size of the clumps. Refer to Tables 
[3] and |4] for details of the clumps identified. 
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Figure 7. The integrated emission maps (contours and grey scale) of HNC and C2H (for the contour levels, see Figure [2} overlaid with 
the corresponding GAUSSCLUMPS clump fits of the 3-D data cube. The ellipses show the orientation and size of the clumps. Refer to Tables 
[5] and |6] for details of the clumps identified. 
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Table 4. Same as Table [s] but for the 20 brightest HCO+ clumps with peak temperatures above Scr level of the Hanning smoothed data 
cube. The clumps are listed in descending peak Tb order. See the online version for a complete list of 186 clumps. 

# RA Dec V Ay Dy Peak Tb dv/dr L Associated clumps in other molecules 



MOUpUl 


lo:zz:Uo.5 


-5U:UD:z8.y 


-49.1 


2.4 


1.7 


1.1 


5.4 


0.6 


23.8 


OoOo, MlNOOl, rllNOlO, OzMOo 


T T * ^\ r\ o 

HCUpUi 


1d:21:U3.2 


rr* or 00 r\ 

-50:35:28.9 


-51.4 


3.8 


1.4 


1.8 


4.9 


0.6 


47.3 


CbOl, Ob02, HJNC02, C2H01 


TT „ n o 


lo:21:25.o 


-50:24:40.9 


-50.1 


2.8 


1.4 


1.5 


4.4 


0.5 


26.1 


ObOo, MlNLyOd, OzMOo 


HC(JpU4 


1d:22:L)9.7 


m rtr' 00 r\ 

-50:06:22.9 


-46.7 


1.7 


1.9 


1.2 


3.0 


0.2 


11.6 


Cb03, Cbl4, HJNCOl, C2H03 


ilL^L^pUO 


iO.Zi.UZ.O 


-ou.oo. iu.y 


-00.0 


O.Z 


n 
u.y 


i .z 


z.y 


n 1 


y.o 


LyoU4, L^ZilUO 


HCOpoe 


16:20:05.8 


-50:57:10.9 


-56.0 


2.0 


1.1 


1.4 


2.9 


0.6 


8.5 


HNC13, C2H20 


HCOpOT 


16:21:14.6 


-50:39:52.9 


-56.0 


3.1 


1.5 


1.9 


2.9 


0.8 


23.6 


CS08, CS17, HNC07, C2H07 


HCOpOS 


16:21:40.3 


-50:23:28.9 


-50.8 


1.9 


1.4 


2.1 


2.6 


0.3 


13.8 


HNC08 


HCOp09 


16:22:11.0 


-50:06:28.9 


-51.5 


2.8 


1.6 


1.2 


2.5 


0.3 


13.6 


CS03, HCOpOl, HNCIO 


HCOplO 


16:20:38.1 


-50:44:16.9 


-53.1 


2.4 


1.2 


1.4 


2.3 


0.4 


8.4 


HNC04, C2H10 


HCOpll 


16:20:44.4 


-50:43:28.9 


-60.5 


1.7 


1.6 


2.9 


2.1 


0.2 


15.1 


HNC14 


HCOpl2 


16:21:36.7 


-50:41:10.9 


-59.8 


2.5 


1.1 


1.4 


2.1 


0.5 


8.0 


HNC09 


HCOpl3 


16:20:39.3 


-50:41:34.9 


-54.6 


2.0 


1.0 


1.4 


2.0 


0.6 


5.1 


CS18, C2H16 


HCOpl4 


16:20:27.3 


-50:41:16.9 


-55.3 


2.3 


1.3 


1.8 


1.9 


0.5 


10.1 


CSIO, HNC16, C2H11 


HCOplS 


16:21:43.5 


-50:28:04.9 


-50.1 


2.0 


1.5 


2.5 


1.8 


0.2 


12.8 


HNC20, C2H17 


HCOpie 


16:21:32.2 


-50:27:04.9 


-52.4 


2.4 


1.1 


1.4 


1.8 


0.3 


6.5 


CS07, CS16, HNC12, C2H02 


HCOplT 


16:20:53.1 


-50:43:58.9 


-54.0 


2.2 


1.2 


2.4 


1.8 


0.6 


10.6 


HNC19 


HCOplS 


16:20:11.5 


-50:53:28.9 


-54.4 


2.6 


1.2 


1.1 


1.8 


0.4 


5.7 


CS05, CS09, HNC17 


HCOpl9 


16:19:48.6 


-51:02:28.9 


-58.9 


2.0 


1.2 


1.5 


1.8 


0.1 


5.6 




HCOp20 


16:21:32.1 


-50:26:22.9 


-49.7 


1.9 


1.2 


2.6 


1.7 


0.3 


10.3 


CS07, CS16, HNC12 



Table 5. Same as Table |3] but for the 20 brightest HNC clumps with peak temperature above 5a level of the Hanning smoothed data 
cube. The clumps are listed in descending T^ order. See the online version for a complete list of 128 clumps. 

# RA Dec V AV Dy Peak Tb dv/dr L Associated clumps in other molecules 



HNCOl 


16:22:09.1 


-50:06:22.9 


-47.1 


3.7 


1 


1 


1.7 


4.4 


0.2 


28.1 


CS03, CS14, HCOpOl, HCOp04 
C2H03 


HNC02 


16:21:02.6 


-50:35:22.9 


-51.7 


3.0 


1 


6 


2.0 


4.2 


0.4 


35.4 


CSOl, CS02, HCOp02, HNC05 
C2H01 


HNC03 


16:21:25.8 


-50:24:46.9 


-50.4 


2.9 


1 


3 


1.6 


3.4 


0.5 


19.8 


CS06, HCOp03, C2H08 


HNC04 


16:20:38.7 


-50:44:22.9 


-52.7 


2.4 


1 


2 


1.4 


2.9 


0.5 


10.9 


HCOplO, C2H10 


HNC05 


16:21:01.3 


-50:35:34.9 


-49.2 


2.1 


1 


4 


1.7 


2.8 


0.2 


13.2 


CS02, HCOp02, HNC02, C2H01 
C2H18 


HNC06 


16:21:02.6 


-50:35:10.9 


-54.6 


2.6 





9 


1.2 


2.8 


0.1 


7.4 


CS04, HCOp02, HCOp05, HNC02 
C2H05 


HNC07 


16:21:15.2 


-50:39:46.9 


-56.1 


3.1 


1 


4 


1.6 


2.4 


1.1 


15.2 


CS08, CS17, HCOp07, C2H07 


HNC08 


16:21:39.1 


-50:23:52.9 


-50.8 


2.3 


1 


4 


2.8 


2.3 


0.1 


19.1 


CS12, HCOp08, C2H15 


HNC09 


16:21:36.7 


-50:41:04.9 


-59.2 


3.1 


1 


2 


1.4 


2.2 


1.1 


10.6 


HCOpl2 


HNCIO 


16:22:09.7 


-50:06:34.9 


-49.7 


2.8 


1 


8 


1.3 


2.2 


0.3 


12.9 


CS03, HCOpOl, HCOp09, C2H03 


HNCll 


16:20:08.9 


-50:53:22.9 


-58.1 


3.3 


1 


9 


2.1 


2.1 


0.6 


25.7 


CS05, CS15, C2H06 


HNC12 


16:21:30.3 


-50:26:52.9 


-52.1 


3.0 


1 


5 


1.5 


2.1 


0.3 


12.6 


CS07, CS16, HCOpl6, HCOp20 
C2H02 


HNC13 


16:20:05.8 


-50:57:04.9 


-56.1 


2.0 


1 


1 


1.0 


2.0 


0.6 


4.3 


HCOp06, C2H20 


HNC14 


16:20:44.4 


-50:43:10.9 


-60.5 


1.7 


1 


4 


2.5 


1.9 


0.1 


10.7 


HCOpll 


HNC15 


16:21:13.9 


-50:33:58.9 


-51.5 


2.0 


1 


1 


1.8 


1.8 


0.4 


7.1 


CS20 


HNC16 


16:20:27.3 


-50:41:28.9 


-56.1 


2.3 


1 


2 


1.7 


1.8 


0.2 


7.9 


CSIO, HCOpM, C2H11 


HNC17 


16:20:11.5 


-50:53:28.9 


-54.6 


2.9 


1 


1 


1.1 


1.8 


0.5 


5.6 


CS05, CS09, HCOpl8 


HNC18 


16:20:47.5 


-50:38:46.9 


-53.1 


2.0 


1 


2 


1.7 


1.8 


0.0 


6.8 


CSll, C2H04 


HNC19 


16:20:51.7 


-50:43:46.9 


-54.2 


2.0 


1 


3 


2.5 


1.7 


0.4 


10.0 


HCOpl7 


HNC20 


16:21:43.5 


-50:28:04.9 


-50.3 


1.9 


1 


2 


2.0 


1.6 


0.1 


6.9 


HCOpl5, C2H17 



the brightness (beam efficiency corrected) of optically thick 
and thin emission lines as follows: 



Tb (thick) 
Tb(thin) 



1 



1 — e~'^thick/-'^ ' 



(1) 



where X is the abu ndance ratio, assuming [CS/C^^S] = 22.5 



( Chin et al.||l996|) and [HCO+/H^^CO +] = 50, based on 



[CO/'^CO] from 



Wilson k Rood 



(1994) 



After solving for r, we next calculate the excitation tem- 



perature Tex from the optically thick emission lines, OS and 
HCO^. From the equation of radiative transfer, Tex can be 
derived from 



Tb = /[J(Te.)- J(Tbg)][l-e-n- 



(2) 



/ is the beam filling factor which we assume to be 1, and 
J(T) = To/(e^o/^ - 1). Once r and Tex are known, the 
column density of the optically thick line can be derived 
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Table 6. Same as Table |3] but for the 20 brightest C2H clumps with peak temperature above 3cr level of the Hanning smoothed data 
cube. The clumps are listed in descending peak Tb order. See the online version for a complete list of 78 clumps. The C2H clumps 
presented here refer to the main hyperfine component only, hence the luminosity is not calculated. 



# 


RA 


Dec 


V 


Ay 




Dy 


Peak Tb 


dv/dr 


Associated clumps in other molecules 


C2H01 


16:21:02.6 


-50:35:34.9 


-50.8 


2.8 


2.2 


1.6 


1.8 


0.5 


CSOl, CS02, HCOp02, HNC02, 




















HNC05, C2H18 


C2H02 


16:21:30.3 


-50:26:52.9 


-52.6 


2.4 


1.5 


1.8 


1.6 


0.2 


CS07, CS16, HCOpl6, HNC12 


C2H03 


16:22:07.8 


-50:06:34.9 


-47.8 


3.2 


1.1 


1.8 


1.6 


0.2 


CS03, CS14, HCOpOl, HCOp04, 




















HNCOl, HNCIO 


C2H04 


16:20:49.4 
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C2H06 
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-58.3 
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1.6 
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C2H07 


16:21:16.5 


-50:39:40.9 


-56.7 


2.3 


1.3 


1.7 


1.1 


0.9 


CS08, CS17, HCOp07, HNC07 


C2H08 


16:21:27.1 


-50:24:40.9 


-50.4 


1.8 


1.6 


1.7 


1.1 


0.5 


CS06, HCOp03, HNC03 


C2H09 


16:19:53.0 


-51:01:34.9 


-47.6 


2.6 


1.0 


1.5 


1.0 


0.4 




C2H10 


16:20:38.7 


-50:44:04.9 


-52.6 


2.0 


1.1 


1.3 


1.0 


0.7 


HCOplO, HNC04 


C2H11 


16:20:26.7 


-50:41:16.9 


-55.5 


2.4 


1.1 


1.8 


1.0 


0.5 


CSIO, HCOpl4, HNC16 
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CS05, CS09, HCOp02, HNC05 


C2H19 


16:21:20.0 


-50:10:04.9 
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16:20:05.8 
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HCOp06, HNC13 



with (cf |Purcell et al.|2006 ) 

N ■■ 



I 



Q{T) / ndv, (3) 



hc^g^A^i 1 - e-^ 

where A^i is the Einstein A coefficient in s""*^, zy is the tran- 
sition frequency in Hz, J T\ydv is the integrated brightness 
temperature in K km s~^, Q(T) is the partition function, 
is the upper state degeneracy and is the upper state 
energy in m~^. At positions, where no C'^'^S or H^'^CO^ was 
detected, we have assumed also that the main isotopologue 
was optically thin and the excitation temperature was set 
to 3 K. Eleven CS and two HCO^ clumps fall into this 
category. The choice of 3 K for the excitation temperature 
was motivated by the median value of Tex from clumps that 



show emission in the main and the rare species (Figures 10 
and 



11) 



Shown in Figure [S] and |9] are histogram plots of CS 
and HCO^ clump column density. The solid line shows the 
column density after correction for optical depth, and the 
dashed line shows the derived column density distribution 
by assuming CS and HCO^ are optically thin with an 
excitation temperature of 20 K. The plots show that after 
correction, the column density increases by an order of 
magnitude. An examination of the clump spectra confirms 
that clumps with column density below 10^"^ cm~^ are 
those with minimal or no optically thin line detections, 
while those with column density of 10^^ cm~^ or above are 
optically thick clumps. 

An excitation temperature of 3 K is low. One possible 
cause could be due to the optically thin lines (C'^^S and 
H^^CO^) being spatially compact, and the assumption 
of having the same beam efficiency as the more common 
isotopologue may not hold. We thus also calculated the 
excitation temperatures of both CS and HCO^ with 
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Figure 8. Histogram of CS clump column density with optical 
depth correction (solid). The dashed line shows the distribution 
obtained by assuming that CS is optically thin and has an exci- 
tation temperature of 20 K. 



extended beam efficiency (?7xb) for CS and HCO^, and 
main beam efficiency (r/nab) for C^^S and H^^CO+. The 
excitation temperatures of CS and HCO^ then increases 
to ^4 K, still low. Another possibility is that the clumps 
are fragmented and unresolvable within the Mopra beam 
(^-^36 arcsec at 100-GHz). The clump finding algorithm 
sees these fragmented clumps as one single clump. Thus, 
the assumption of the beam filling factor (/) being unity 
does not hold, despite the fact that all identified clumps 
are larger than 1.5 beam widths. Various studies show 
that CS clumps are typically fragmented, with a volume 
filling factor as low as 0.2 (e.g. |Stutzki & Guesten 1990 



puvela|1998| ) . To investigate the effect of beam filling factor 
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Figure 9. Histogram of HCO"^ clump column density with opti- 
cal depth correction (solid). The dashed line shows the distribu- 
tion obtained by assuming that HCO"^ is optically thin and has 
an excitation temperature of 20 K. 
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Figure 10. Histogram of CS clump excitation temperature with 
a beam filling factor of 0.2 (solid line). The dashed line shows the 
distribution assuming a unity beam filling factor. 
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Figure 11. Histogram of HCO"^ clump excitation temperature 
with a beam filling factor of 0.2 (solid line). The dashed line shows 
the distribution assuming a unity beam filling factor. 
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Figure 12. Plot of variances of the principal components. The 
inset presents a zoomed version of the plot, showing the second 
and higher order components in more detail. 



on the excitation temperature, we derived the excitation 
temperatures of CS and HCO^ clumps with various beam 
filling factors. We found that with / = 0.2, approximately 
74 per cent of the CS clumps, and 56 per cent of the HCO^ 
clumps have Tex between 6 to 20 K. Histograms of CS and 
HCO^ excitation temperatures are shown in Figure 



10 



and 

|ll| to illustrate the effect of lowering beam filling factor. The 
mean excitation temperature increases greatly by lowering 
the beam filling factor, suggesting the clumps may indeed 
fragmented and not resolvable with MOPRA beam. This 
would provide an explanation for the lack of correlation be- 



tween clump size and line width (see Section 6.1 for details) 



5 PRINCIPAL COMPONENT ANALYSIS 

An alternative way to study the distribution of the 
molecules is by performing principal component analysis 
(PCA) on the integrated emission maps of the region. PCA 
is a multivariate data analysis technique, its purpose being 
to reduce the dimensionality of a data set. Mathematically 
the derivation of principal components involves finding the 
eigenvalues and eigenvectors of a covariance or correlation 



matrix (s ee e.g. |jolliffe| |2002'; "Hey er Schloerb] [1997 
[Ungerechts et al.||1997 and references therein) 



The purpose of performing PCA on the integrated 
emission maps is to characterise differences in molecular 
distribution. The G333 cloud consists of star forming sites 
at different stages of evolution, from cold dense starless 
cores to Hll regions. Hence we expect to see chemical 
differences across the cloud. The data set covers varieties 
of chemical probes. Representative ones are the outflow 
tracer HCO^, the large scale gas tracer ^^CO, the dense gas 
tracer CS, and the cold dense gas tracer N2H^. Interesting 
aspects are chemical differences and a possible temperature 
dependance of HCN to HNC abundance ratio. In this 
section, we will discuss the result of a PCA decomposition 
on the whole G333 molecular cloud; the PCA of individual 
regions of interest will be discussed in later sections. 

We have chosen eight molecule lines with high signal- 
to-noise ratio for this analysis: ^^CO, C^^O, CS, HCO+, 
HCN, HNC, N2H+ and C2H (Figures [l] and |2|. Among 



16 A^. Lo et al 




16^22°'30^ ZZ'^Q^ 30^ Sl^'O^ 30^ 201^0^ 30^ le^'SS^'SO^ SS^'O^ 30^ 2l"^0^ 30^ SO^'O^ 30^ 

RA (J2000) RA (J2000) 



Figure 14. Constructed images of the first four principal components of the GMC. The contour levels are multiples of 20 per cent of the 
peak. The dashed contours represent an anti-correlation, while the solid contours represent a positive correlation in each of the principal 
components. The grey scale in each of the images shows the first principal component, overlaid with contours of the other principal 
components for ease of comparison. The lower left circles of each images indicate the beam size. The arrow indicates an example of 
anti-correlation between HCO+ and N2H+ as discussed in Section [6.2.2] 
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Table 8. The correlation matrix of the input molecular data set for principal component analysis. 

CS HCO+ HNC C^^O C^ii HCN N2H+ ^^CO 



cs 


1.00 












HCO+ 


0.85 


1.00 










HNC 


0.90 


0.92 


1.00 








CiSQ 


0.77 


0.80 


0.81 


1.00 






C2H 


0.60 


0.70 


0.67 


0.61 


1.00 




HCN 


0.81 


0.94 


0.90 


0.75 


0.70 


1.00 


N2H+ 


0.80 


0.76 


0.82 


0.73 


0.56 


0.71 


13CO 


0.77 


0.82 


0.81 


0.88 


0.63 


0.79 



Table 9. The eigenvectors and eigenvalues of the principal components (PC) derived from the correlation matrix listed in Table [8^ 





Percentage of variance 


CS 


HCO+ 


HNC 


CiSQ 


C2H 


HCN 


N2H+ 


13CO 


PC 1 


80.0 


0.36 


0.38 


0.38 


0.35 


0.30 


0.37 


0.34 


0.35 


PC 2 


6.3 


-0.25 


0.05 


-0.10 


-0.17 


0.86 


0.15 


-0.36 


-0.06 


PC 3 


4.7 


-0.25 


-0.11 


-0.22 


0.58 


-0.07 


-0.16 


-0.38 


0.61 


PC 4 


3.8 


-0.05 


-0.33 


-0.17 


0.21 


0.38 


-0.50 


0.64 


-0.03 


PC 5 


2.1 


0.80 


-0.26 


0.07 


-0.05 


0.15 


-0.34 


-0.39 


0.04 


PC 6 


1.4 





0.03 


0.13 


0.68 


0.01 


0.02 


-0.18 


-0.70 


PC 7 


0.8 


-0.25 


-0.63 


0.72 


-0.03 


0.01 


0.13 


-0.05 


0.09 


PC 8 


0.7 


0.22 


-0.52 


-0.48 


0.10 


-0.01 


0.66 


0.11 


-0.03 



these, C2H has the lowest signal-to-noise level and ^^CO 
has the highest. The data sets were imported into idl 
(a data visualisation and analysis platform) and each of 
their standard deviations were obtained. We decided to 
derive the principal components of the correlation matrix, 
rather than the covariance matrix. This is to avoid the 
principal components being dominated by a single variable 
(Jolliffe 2002 Chapter 3.3). An alternative approach 
involves normalising the data set first, and then deriving 
the principal components from the covariance matrix. We 
in fact derived the principal components using these two 
different approaches and obtained essentially the same 
result. 

Listed in Table [S] is the correlation matrix of the input 
molecules, which describes how well the molecules correlate 
with each other. The correlation matrix shows there are 
four pairs of molecules which have correlation coefficients 
equal to or above 0.9. These are HCO+ and HCN, HCO+ 
and HNC, HCN and HNC, HNC and CS. This is consistent 
with the integrated intensity maps shown in Figure [l] and 
U where CS, HCO+, HCN and HNC have similar large 
scale spatial distribution. However CS does show small 
differences in distribution, as discussed in Section |3.1| 
Among the eight molecules, C2H and N2H^ have the lowest 
correlation coefficient. After forming the correlation matrix 
of the data set, we derived the eigenvalues and eigenvectors 
of the correlation matrix (Table [9|. The eigenvalues indicate 
that the first principal component accounts for 80 per cent 
of total variation, while the first three principal components 
together account for over 90 per cent of total variation 
in the data set (see Table [9] and Figure 12). Figure 12 



also shows that the first four principal components contain 
features above the noise level, while the variations contained 
in higher components are insignificant. 

From the eigenvectors we constructed images of 
principal components (PC), by projecting the data set 
onto each of the eigenvectors. Shown in Figure [13] are 
plots of eigenvectors of each molecule in the first four 
principal components. A (negative) positive value indicates 



the molecule is (anti-) correlated with others. The larger 
the value, the stronger the correlation. From the PCI 
axis of Figure |13^ it becomes clear that all of the eight 
molecules are positively correlated with each other. This 
can be visualised in Figure |14^ , which is the whole data 
set projected onto the first principal component only. This 
also resembles an 'ideal' molecular distribution of the giant 
molecular cloud - the average intensity distribution of the 
species - and can be interpreted as the eight molecules being 
positively correlated on large scales. From the eigenvector 
plot of the second principal component (PC 2 axis of Figure 



13 1), C2H stands out from other molecules; we suspect this 



is caused by scanning stripes resulting from the on-the-fiy 
mapping procedure. This can be seen in the image of the 



second principal component (Figure 14 3); note the vertical 



and horizontal stripes of solid contours which resemble 
the scanning patterns presented in the C2H integrated 
emission map (Figure [2]l). In fact, the scanning stripes 
are present in all data taken simultaneously as listed in 
Table [1] (the July, 2006 observation season). However, due 
to the lower signal-to- noise level of the C2H emission, these 
scanning patterns become dominant in this map. We have 
explored the possibility of removing scanning artefacts by 
reconstructing the images without the second principal 
component. This indeed reduced the level of the artifacts, 
but it also subtracts a component of real emission from the 
image. Hence it provides a qualitative improvement, but 
the corrected image should not be used for quantitative 
studies. 

From the eigenvector plot of the third principal 
component (Figure 13 d) we can see C^^O and ^^CO are 
anti-correlated with other molecules. Spatially this anti- 
correlation appears as compact negative contours against 
a diffuse positive region in the third principal component 



image (Figure 14:). This could be due to the differences 
between high and low density tracers discussed in Section 
|6.2.1| with the CO isotopologue distribution being more 
extended. 

The fourth principal component eigenvector plot 
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6 DISCUSSION 
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In this section we discuss the results from this multi- 
molecular line survey and the analysis techniques applied. 
We also focus on particularly interesting regions in the 
G333 giant molecular cloud. 
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Figure 13. Plots of eigenvectors of the first three principal com- 
ponents of the data set. Each of the eigenvectors represent the 
component of that molecule in the relevant principal component. 



(Figure 
HCO+. 



13 



3:) shows an anti-correlation between N2H^ and 
The anti-correlation also appears clearly in the 
fourth principal component image (Figure 141), notably 
across IRAS16172— 5028 (shown by the arrow), possibly 
pointing to chemical differences across the object. Further 
discussion on this will be presented in Section [6. 2. 2 1 

The principal component analysis provides a way to 
quantify the spatial distribution of different molecules (with 
the constructed principal component images), and in plac- 
ing different molecules into groups (with the eigenvectors). 
This results in minimising the number of different molecules 
required to characterise the basic properties of the cloud: 
CO isotopologues, HCO^, N2H+, and C2H. Implications 



6.1 Size-line width, luminosity-size, and 
luminosity-line width correlation 

To examine the correlation among the clump properties of 
CS, HCO^ and HNC, we have compared luminosities of the 
molecular line emission, line widths and radii in Figure [15] 
The C2H clumps presented here were obtained analysing 
the main hyperfine component only. Hence the entire line 
luminosity has not been calculated. The luminosity was cal- 
culated with 



{d[pc]f{ 



180 X 3600^ 



•I 



Tdv, 



(4) 



where d is the distance to the CMC (3.6 kpc), rx and ry 
are the radii (arcsecond) of the two principal axes and 
f Tdv (K km s~^) is the sum of emission at the maximum 
position. We have discarded clumps with line widths below 
0.7 km s~^, as this is less than twice the effective spectral 
resolution of 0.3 km s~^ after hanning smoothing and 
binning the data. The slope of the ordinary least square 
(OLS) bisector ( [Feigelson L Babu|1992t line of best fit (7), 
estimates the structural relationship between variables X 
and y, without assuming whether Y depends on X or vice 
versa. Values inside the curly brackets are the slopes of 
OLS (X|y), which minimises the residuals in X, and OLS 
(y|X) which minimises the residuals in Y. The number of 
clumps (N) are shown on the top left corners of correlation 
plots. 



Larson ( 1981 ) empirically showed, utilising data from 



several molecular cloud surveys, that the line width AV 
appears to scale with cloud size r, AV oc r^. The com- 



mon value quoted for 7 is ^ 0.4 ± 0.1 (see [Mac Low 
Klessen 2004, and references therein). However, significant 
deviations have been reported (e.g. |Loren| |1989| [Caselli 



Myers] |1995| [Plume, Jaffe, Evans II, Martin- Pintado 



& Gomez- Gonzalez 



1997) 



Consistent with our analysis 
of '^CO and C'^O^ihir^mission (WLB2008), we find no 
correlation between line width and radius for CS, HCO^ 
and HNC clumps. The line width-size (note radius is given, 
not diameter) plots (Figure 15i-c) show a large scatter, 
hence the fitted slopes are poorly constrained. Note the 
large differences between OLS (-^|1^) and (1^|-^) slopes. 
Other line width-size relation studies (e.g. [Schneider 
Brooks 2004| and references therein) also found no line 



width-size relation, and suggested its possible dependance 
on clump identification procedure. 

There is a clear correlation between luminosity and 
radius among CS, HCO^ and HNC clumps. The fitted 
OLS bisector slopes 7 (3.8 to 4.6) are slightly steeper 
than found for ^^CO and C^^O (- 3.0). Since luminosity 
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Figure 15. Plots of relationship between physical parameters for CS, HCO+ and HNC clumps: (a) to (c) radius and line width; (d) to 
(f) luminosity and radius; (g) to (i) luminosity and line width. 7 indicates the slope of the ordinary least square (OLS) bisector fitted 
solid line, the values in curly brackets are the slope of the OLS and (y|X). N is the number of clumps. Note that clumps with 

line widths smaller than 0.7 kms ~^ are discarded. 



implicitly depends on r^, one would expect 7 ~ 3, assuming 
a constant column density across the clumps. However, no 
statistical significance can be assigned to the value of slope, 
due to the inherent bias in clump finding algorithms, that 



tend to find larger clumps in denser regions ( [Schneider 
Brooks|2004 



Similarly, luminosity and line width exhibits a strong 
correlation; this is not surprising as luminosity is calculated 
from the clump sizes and line widths. We do note that the 
^^CO and C^^O clumps from previous work (WLB2008) 
were found using a different clump- finding algorithm, 
CPROPS. 



6.2.1 High and low density tracers 

The eigenvector plot of the third principal component shows 
that ^^CO and C^^O are anti-correlated with the other 



molecules (Figure 13 d); this can be visualised spatially from 
the projected principal component image as shown in Figure 
14:. In this third PC image the compact dashed (negative) 
contours contrast spatially with the diffuse solid (positive) 
contours. This is due to C^^O and ^^CO (J = 1 ^ 0) 
being low density tracers, and other molecules presented 
here being high density tracers (see Table [l] column 5). 
^"^CO and C^^O are present in both low and high density 
regions, whereas the other molecules are present in dense 
gas only. Hence in diffuse gas ^^CO and C^^O are in 'excess' 
compared to the high density tracers. 



6.2 Physical and chemical properties from PCA 

In Section [5] we presented the results of principal component 
analysis, finding correlations or anti-correlations between 
molecular species. In the following sections, we discuss these 
in terms of the physical and chemical properties of the GMC. 



6.2.2 The anti- correlation of N2H^ and HCO^ 

The eigenvector plot for the fourth principal component 
(Figure 
HCO+ 



13 



:) indicates the two ionic species N2H^ and 
are anti-correlated. According to models of dense 
clouds the major formation pathway o f N2H^ is via ion- 
molecule reactions between N2 and HJ (Nejad, Williams & 
Charnley[l9 9Q). It is mainly destroyed by reactions with CO 



^ Tafalla, Myers, Caselli &: Walmsley|2004 ) and by recombi- 
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nation with electrons in hot regions ( Sternberg & Dalgarno 



1995 ). The following chemical equations summarise the men- 



tioned reactions, 



N2 + m 



N2H+ + H2 , 



N2H+ + CO — > HCO+ + H2 , 



(5) 
(6) 
(7) 

HCO^ can be formed by ion-molecule reaction between 
N2 H+ and CO (Equation § 



or between CO and HJ (jStern- 



berg k, Dalgarno 1 995) as shown in Equation [s] below, 

CO + H+ — > HCO+ + H2 . (8) 

In cold dense cores, CO, the major destroyer of N2H^, will 
be depleted, so that the N2H^ abundance increases. With 
the lack of CO molecules, HCO^ loses both of its major 
pathways (CO reacting with N2H^ or HJ), so the abun- 
dance of HCO^ decreases. On the other hand, in warm re- 
gions where CO is not depleted, it reacts with N2H^ and 
produces HCO^ causing a decrease in N2H^ and increase 
in HCO^. CO reactions with further increase the HCO^ 
abundance. 

Prominent sites to study this phenomenon are located 
in the vincinity of Hii regions, i.e. in photodissociation or 
photon dominated regions (PDRs), where at the cloud sur- 
face a high HCO^ and a low N2H+ abundance is expected. 
In contrast, in the inner part of the cloud, where material 
is collapsing to form a dense core, CO is depleted causing 
an increase in N2H^ and a drop in HCO^ abundance. Thus 
HCO^ and N2H^ should be anti-correlated. 



6.2.2.1 IRAS16172-5028 The Hii region associated 
with IRAS16172— 5028 provides a good candidate to test 
whether N2H^ and HCO^ are anti-correlated. As suggested 
by the compact intense N2H^ emission, this source has a 
cold high density component, but is immersed in a strong 
FDR according to the bright 8.0-/im emission. To exclude 
the effects of large-scale spatial contamination, we applied 
PCA to the integrated emission map of this region. Listed 
in Table [To] are the eigenvalues and eigenvectors for this 
analysis; the eigenvectors of the second, third and fourth 
principal components are also plotted in Figure |16| where 
the anti-correlation of HCO^ and N2H^ is clearly seen. 
In the eigenvector plot of the second and third principal 
components (Figure 16 1), C2H has the strongest positive 
correlation among the molecules. This is due to its spa- 
tially wide-spread emission towards the south-eastern re- 

N2H+. 
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shows the 



gion, compared to, for example, N2H^. Figure 
GLIMPSE 8.0-/im image (grey scale) overlaid with contours 
of the second principal component of the Hii region associ- 
ated with IRAS16172— 5028, demonstrating that there are 
two regions with molecular emission that are anti-correlated. 
An examination of the eigenvectors shows that C2H, HCO^ 
(thick contours) and N2H+ (thin contours) contribute sig- 
nificantly to the variance in the second component and are 
anti-correlated, as expected for the hypothesis above. Note 
that the 8.0-/im emission (PAHs and warm gas) coincides 
well with the thick contours (HCO^). 

The second and third principal components of this re- 
gion also show that ^^CO and C^^O are correlated with 





1.0 




C2H ■ 




0.5 


■ 












Q_ 


0.0 




r\HCO+ 




-0.5 




cs 




-1.0 




.. 1 .... 1 .... 1 ... . 




-1 


.0 


-0.5 0.0 0.5 1. 




1.0 




PC 2 


(b) 




T . 1 . 1 . . 1 . 1 . . 1 . 1 . 1— 

CS 




0.5 




\ N2H+ 




D. 


0.0 




\ / ^^C2H 




-0.5 




HCN " 




-1.0 







-1.0 -0.5 0.0 0.5 1.0 
PC 3 



Figure 16. The eigenvectors of principal components two, 
three and four of the gas near the Hll region associated with 
IRAS16172-5028. 
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Figure 17. An image of the second principal component 
(contours) of the gas near the Hll region associated with 
IRAS16172-5028 overlaid on the GLIMPSE 8.0-/im image (grey 
scale). The thin contours show regions with negative correlation 
to this component (mainly N2H"^) and the thick contours are 
regions with positive correlation (mainly HCO"^). 
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Table 10. The eigenvectors and eigenvalues of the principal components (PC) of the gas near the Hll region associated with 
IRAS16172-5028. 
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Table 11. The eigenvectors and eigenvalues of the principal components (PC) of the gas near the Hll region associated with 
IRAS16177-5018. 





Percentage of variance 
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Figure 18. The eigenvectors of principal components two, 
three and four of the gas near the Hll region associated with 
IRAS16177-5018. 



N2H^ (Figure 16 1), which contradicts the idea suggesting 
that CO is the main destroyer of N2H^. However, we be- 
lieve this can be explained by considering the different den- 
sity conditions that these two species are tracing. CO is a 
low density tracer having a critical density of 10^ cm"*^, 
while N2H^ is a high density tracer with a critical den- 
sity of ~ 10^ cm~^. Therefore CO is tracing both low and 
high density gas along the line of sight. Since maps of inte- 
grated emission were used for the PCA, the ^^CO and C^^O 
may arise from a molecular envelope almost devoid of N2H^, 
while the bulk of the N2H+ emission is arising from the core. 



6.2.2.2 IRAS16177-5018 Another example is the 
molecular emission associated with IRAS16177— 5018 as 
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Figure 19. An image of the third principal component (contours) 
of the gas near the Hii region associated with IRAS16177— 5018 
overlaid on the GLIMPSE 8.0-)Um image (grey scale). The thin 
contours showing regions with negative correlations (mainly 
N2H"^) and the thick contours are regions with positive corre- 
lations (mainly ^^CO and C^^O). 



shown in Figures ^] and |20| The eigenvalues and eigen- 
vectors of the principal components are listed in Table [ll] 
while the eigenvectors of the third and fourth principal 
components are plotted in Figure |18| Similar to the other 
source IRAS16172— 5028, this region also shows an anti- 
correlation of N2H^ and HCO^ in the fourth principal 



component (Figure 20). On the other hand, shown in the 
third principal component (Figure 19 ) is an anti-correlation 
of N2H^ and CO isotopologues. 

The PCA also indicates that HCN and HCO+ are 
correlated, as shown in the eigenvector plots of the entire 



22 A^. Lo et al 



-50°22' 




1500 



1000 



500 



16^21^50^ 40^ 30^ 20^ 10^ 
Right Ascension (J2000) 



Figure 20. Similar to Figure ^9] but the fourth principal com- 
ponent (contours) overlaid on the GLIMPSE 8.0-/im image (grey 
scale). The thin contours show regions with negative correlation 
to this component (mainly HCO"^) and the thick contours are 
regions with positive correlation (mainly N2H"^). 



region (Figure 13). This correlation has also been suggested 



by other studies, such as by Turner & Thaddeus (1977) 
for Orion-KL. IRAS16172-5028 and IRAS16177-5018 are 
both known to be associated with an Hii region, contain 
cold molecular gas of high density (bright compact N2H^ 
emission), and have a low photon exposure. As expected by 
chemical models, both sources clearly show that N2H^ and 
HCO^ are anti-correlated. 



6.3 The infrared ring 

Another note worthy feature in G333 is a ring of mid- 
infrared emission as shown in Figure [21] the northern por- 
tion of which is aligned with an arc of N2H^ emission (see 
the top arrow on Figure |2|. This ring feature is visible in 



the GLIMPSE 5.8 and 8.0-/im (Figure [21j top right and 
bottom left panels) images, but is marginally detectable 
(Figure 



in 3.6- 
GLIMPSE 5 
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and > 



top left panel) and 4.5-/im images. The 
.0-/im channels of the Spitzer IRAC in- 
strument are dominated by PAH emission excited by nearby 
ultraviolet sources (Reach et al. 2QQ6| ). Lying inside the 



emission ring but offset towards the edge, is bright MIPS- 
GAL ( Carey et al.|2QQ6 ) 24- /xm emission (Figure [2l] bottom 
right panel) , spatially coincident with the MSX point source 
G333 . 5114-00. 2798 (marked with a white circle; Egan et al. 
2003|) and an 5' source IRAS16182-5005 (marked with 



triangle). The absence of the emission ring at 24- /xm further 
suggests that the emission originates from PAHs, and that 
there is a lack of thermal dust emission. 

Comparing the MIPSGAL 24-/im emission with the 
843-MHz radio continuum (Figure 211 contours) from the 
Molonglo Galactic Plane Survey (MGPS, [Green et al.|1999t , 
it is clear that the radio continuum is associated with the 
24- /im emission. A search of the SIMBAD astronomical 



database and VizieR catalogue servic^did not lead to the 
identification of any supernova or X-ray source, suggesting 
there is an Hii region inside the 8-fim emission ring. Shown 
in Figure 21 j is the 5.6-/im emission ring (grey scale) over- 
laid with 1.2-mm dust continuum contours (Mookerjea et al 



2004); the 1.2-mm dust continuum lies next to the bright 8- 



fim rim and is coincident with a dark filament (appearing 
white in the grey scale image). According to 'Mookerjea e t ah] 
(|2004|), this source (MMS16) has a density of - 2x10^ cm"''. 



The 1.2-mm dust continuum em ission agrees well with the 
arc of N2H^ emission (Figure 21 :) described above. There is 
no detectable dust continuum from the south-western region 
of the ring. Molecular emission is detected from this part 
but it is weak and diffuse, indicating a relative low density 
region compared to the dust. Given the completeness and 
symmetry of the ring, one might expect the Hll region to 
lie near the centre. Since it does not, we suggest the possi- 
bility that this is due to density differences in the gas sur- 
rounding the infrared ring. In general Hii regions are density 
bounded, and the pressurised Hii gas breaks out of the cloud 
into lower-density gas, creating a champagne flow (Stahler 



& Palla 2005). The inhomogeneous density (from our molec- 
ular line data) allows the ionised gas to spread out further 
towards the south-west compared to the north-east rim, cre- 
ating a near circular 8-/im emission structure with the driv- 
ing source lying on the edge. However the extreme symmetry 
also suggests the ring could be a pre-existing structure that 
is currently being illuminated by the Hii region. In either 
way this infrared emission ring must be related to the Hii 
region. 

Another notable feature is the spatial coincidence of 
mm dust continuum, N2H^ emission and the infrared dark 
filament (which appears as white in the grey scale image) to 
the north of the ring (Figure 21 3 and c). As CO is frozen 



on to grains in dense cold gas (see e.g. Kramer et al.[ 1999 



Bacmann et al. 2002), causing it to be depleted, and the 



ion-molecule formation reaction for N2H^ proceeds regard- 
less of temperature, N2H^ is a good candidate for tracing 
cold dense gas where CO is depleted. With mm dust contin- 
uum tracing cold dust along with N2H^ tracing cold dense 
gas, it is not surprising that they coincide with the infrared 
dark filament. 



7 SUMMARY 

In this paper we have presented data from a multi-molecular 
line survey of the southern star forming region, the G333 
giant molecular cloud complex. For this survey, we have ex- 
ploited the Mopra MMIC receiver and the 8-GHz bandwidth 
UNSW-Mopra Spectrometer, resulting in over twenty multi- 
molecular transition maps, with velocity resolution of — 0.1 
km s~^. We have presented total intensity maps of molecules 
with bright emission and have also discussed the velocity 
structure of the G333 molecular cloud. To further charac- 
terise the physical and chemical properties, we have carried 
out common analysis techniques such as GAUSSCLUMPS to 
obtain distributions of CS, HCO^, HNC and C2H emission. 
We have also performed principal component analysis on 
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Figu re 21. (a) GLIMPSE 3.6-/im image; (b) GLIMPSE 5.8-/im image overlaid with 1.2-mm dust continuum contours ( [Mookerjea et al.] 
|2Q04| >; (c) GLIMPSE 8.0 -jim images overlaid with N2H"^ integrated emission contours and (d) MIPSGAL image overlaid with 

MGPS 843-MHz radio continuum contours of the infrared ring discussed in Section [6^ The grey scale is stretched to - 20 MJy sr~^ 
for the GLIMPSE 3.6-/im, - 120 MJy sr-^ for the 5.8-/im image, - 300 MJy sr-Mor the 8-/j,m image and - 800 MJy sr-^ for the 
MIPSGAL image. The contours are at multiples of 20 per cent level of each of the emission peaks. The hatched circles represent the beam 
size of the respective contour maps. The small circle indicates the position of the MSX point source G333. 5114-00. 2798, the triangle 
marks the position of the IRAS source, IRAS16182— 5005 (resolution of ~ 4 arcmin) and the cross shows the position of a far-infrared 
source ( [Karnik et al.|20Ql[ resolution of ~ 1 arcmin) which peaks at 100-/im. Note the ring feature present at 5.8 and 8.0-/im, which is 
dominated by PAHs, the bright 24-^m emission which lies inside this ring. Both the 1.2-mm dust continuum and the integrated N2H+ 
emission coincide with the infrared dark filaments, which appear white in the 5.8 and 8.0-/im grey scale images. 



the data set, to visualise and parameterise the differences 
between the spatial distribution of molecules. In this work, 
we have found: 



(i) Differences in spatial and velocity distribution among 
different molecules. We found that the spatial distribution 
of CS, HCO^, HON and HNC are similar on large scales, 
while N2H^ seems to trace preferentially the very densest 
regions. C2H is only detected close to bright infrared emis- 
sion regions. The detected molecules all have similar velocity 
distributions. 

(ii) The velocity gradient across the GMC complex noted 



in ^^00 (BWC2006) and C^^O (WLB2008) is also present 
in CS, HCO+ and HNC. 

(iii) CS, HCO^ and HNC emission maps were decom- 
posed with GAUSSCLUMPS in three dimensions. We found 
no correlation between clump radius and line width, but a 
clear correlation between luminosity and radius. Account- 
ing for saturation effects in the CS ( J = 2 ^ 1) and HCO^ 
( J = 1 ^ 0) lines toward clumps, we obtain column densi- 
ties of - 10^^ to - 10^^ cm-^ 

(iv) An alternative approach used to characterise this 
data set was principal component analysis (PCA). PCA 
separates molecules into low (^^CO and C^^O) and high 
(the rest) density tracers, identifies anti-correlations be- 
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tween HCO^ and N2H^, correlations between HCN and 
HCO^, and helps to explore scanning patterns of the 'on- 
the-fly' mapping. 

(v) A noteworthy 'ring-like' structure of the GMC is 
present in the GLIMPSE 8-/im image, associated with an 
Hll region as suggested by radio continuum emission inside 
the ring. The molecular line data (especially N2H^) shows 
gas being swept up and compressed by the ring. 
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